Structure and transport in multi-orbital Kondo systems 
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Abstract 

We consider Kondo impurity systems with multiple local orbitals, such as rare earth ions in a metallic host or 
multi-level quantum dots coupled to metallic leads. It is shown that the multiplet structure of the local orbitals 
leads to multiple Kondo peaks above the Fermi energy Ef, and to "shadow" peaks below Ef- We use a slave boson 
mean field (MF) theory, which recovers the strong coupling Fermi liquid fixed point, to calculate the Kondo peak 
positions, widths, and heights analytically at T = 0, and NCA calculations to fit the temperature dependence of 
high-resolution photoemission spectra of Ce compounds. In addition, an approximate conductance quantization for 
transport through multi-level quantum dots or single-atom transistors in the Kondo regime due to a generalized 
Friedel sum rule is demonstrated. 
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1. Introduction 

In realistic Kondo impurity systems [I] the local 
magnetic moment is often distributued over several lo- 
cal orbitals. For example, in rare earth ions, embedded 
in a metallic host, the local moment is carried by the 
seven 4f orbitals whose degeneracy is lifted by spin- 
orbit (SO) and crystal field (CF) splitting. Small quan- 
tum dots (Qdot) containing few electrons are in gen- 
eral comprised of several single-particle levels and can 
be seen as artificial atoms. Such systems, coupled to 
metallic leads, can be tuned to the Kondo regime as 
well, if the Coulomb blockade is large enough to fix the 
electron number on the Qdot to an odd number and, 
hence, the effective dot spin to 1/2. In rare-earth sys- 
tems Kondo resonances have recently been observed 
by high-resolution photoemission spectroscopy [2,3], 
while the advances in nanotechnology have made it 
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possible to perform direct tunneling spectroscopy of 
the Kondo resonance in Qdots [4,5,6]. 

In this contribution we show that each of the 
low-lying single-particle orbitals generates a many- 
particle Kondo resonance near the Fermi energy Ef, 
and we calculate their positions, widths and heights 
at temperature T = 0. The results of a Non-Crossing 
Approximation (NCA) calculation at finite T are 
compared to high-resolution photoemission spectra 
of the rare-earth Kondo system CeSi2- Moreover, 
an approximate conductance quantization in multi- 
orbital Kondo Qdots is predicted originating from a 
generalized Friedel sum rule [7] . 

The systems discussed above can be described by the 
multi-orbital Anderson impurity Hamiltonian 
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where Hm, 
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describes the conduc- 
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tion electron band, and c( m<T , creates an electron with 
spin a in the low-lying single-particle level Edm < 0, 
m = 0, ...,M — 1, with Edo the (noninteracting) 
ground state. Electrons in any of the local orbitals 
have a Coulomb repulsion U and are coupled to 
the conduction band via the hybridization matrix 
elements V mp - For later use we indroduce the effec- 
tive couplings T mm / = ■n^2 p V mp V* rn ,A p<T {Q), with 
A ptT (u}) the conduction electron spectral function and 
N(oj) = A pa (u)) the density of states per spin. In 
the following we will assume the system to be in the 
Kondo regime, T mm /\edm\ *C 1, U — » oo, and for the 
expansions leading to Eqs. (7), (11), and (13) below, 

A mm /(Sdm - £do) < I- 



2. Spectral features 

In 2nd order perturbation theory, the spin scatter- 
ing amplitude of the local spin shows log divergences at 
the energies (e dm - e d0 ) < |edo| (*) and -(edm - £do) 
(**). These instabilities can be understood as due to 
spin fluctuations involving a transition of the local elec- 
tron from the ground state Edo to an excited state Edm 
(m > 0) (*) or vice versa (**) [3,8]. The positions of 
these resonances at T = can be estimated analyt- 
ically from a slave boson MF theory [9] for Eq. (1), 
since this approach recovers the strong coupling Fermi 
liquid fixed point [1] of the model. The numerical eval- 
uations are done within the NCA [3], which correctly 
describes the life-times of the excited resonances and 
the temperature dependence for T ~ Tk [10]. 

In order to implement the effective restriction of no 
double occupancy of all orbitals, enforced by the strong 
Coulomb repulsion U, the local electron operator is 
represented in terms of pseudofermions f ma and slave 
bosons b as d ma = f m<y b, with the operator constraint 
Q := J2ma fma fma + b 1 b = 1. In a functional integral 
approach, the latter is imposed by the functional Kro- 
necker delta 
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Introducing mean field and fluctuating parts r, a for 
the Bose field, b = r + a, ¥ = r + a' and A , A for the 
constraint field, A = —i\ + A, Eq. (1) with Eq. (2) 
takes the form [1] 

^ ^ ^dmfmafma H~ ^ ^ [Vmp/mg^pg 

+ h.c] 
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Fig. 1. Diagrammatic representation of a) the fluctuating Bose 
sclfcncrgy and b) the local physical electron Green's function. 
The solid circle represents the modulus squared of the boson 
MF. 
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Here we have defined renormalized local levels idm 
such that Ao = ido — £do and idm = £dm + Ao and effec- 



tive hybridizations V„ 
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mean field solutions r, ido are then obtained by min- 
imizing the free energy, which is equivalent to setting 
the terms linear in the fluctuating fields a, a\ A equal 
to zero, 



r 2 + (a t o)+^(/+ lCT / mCT )-l = 

ma 
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The expectation values in Eqs. (5), (6) are calculated 
from the Fourier components of the pseudofermion, 
the mixed fermion-conduction electron, and the 
fluctuating boson Green's functions Gf mm i „(t) = 

-l(Tfma(t)fL, a (0)),G c f pma(t) = -i (f C pa (t) fL (0)> , 

G a (t) = -i(fa (t)a^(0)}. Calculating the / selfener- 
gies for a flat band and keeping only the leading terms 
in the effective hybridizations, they are [8,11], 

Gf mma(i^) = [iw — idm — if mm] (7) 
G c f pma(iu) = Vmp G c pa («*>) Gf mma(i^) (8) 

G a (iv) = [if - (£do - £do) - T, a (iv)] 1 . (9) 

As seen from Eq. (14) below, the low-energy peak in 
each of the pseudofermion spectral functions G / mm is 
carried over to a quasiparticle peak in the physical im- 
purity electron spectral function, i.e. to a Kondo res- 
onance. Hence, the width r mm of G/ mm defines the 
Kondo temperature of the mth Kondo peak, r mm — 
Tj(m- It is to be determined from the MF solution be- 
low. The excitied resonances have an additional width 
due to inelastic relaxation [10]. The selfenergy of the 
fluctuating Bose field appearing in Eq. (5) is given by 
(Fig. 1 a)), 
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Fig. 2. Local electron spectrum of the Hamiltonian Eq. (1), 
NCA calculation, M = 3. e d0 ,i,2 = -0.95-D, -0.6D, -0.4D, 
r 00 = Tn = r 22 = 0.2D; T/co = 3 ■ 10 _3 i3; T = 5 ■ 10~ 5 D; 
D— half band width. 

It is worth noting that Eq. (5) just represents the 
Friedel sum rule, when the expectation value (a* a) of 
the fluctuating boson density is included [8,11]. There- 
fore, it describes the correct shift of the Kondo reso- 
nances w.r.t. Ef due to potential scattering. We now 
use the relations i d0 < T K o and i dm > T Km , m = 
1, . . . , M — 1, where the former will be justified below 
by the MF solution in the Kondo regime, and the latter 
is true for sufficiently large splitting of the local levels. 
The positions of the Kondo peaks and the renormalized 
Tk are obtained from Eqs. (6), (5), respectively, as 
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Finally, we obtain the physical electron Green's func- 
tion by evaluating the diagrams shown in Fig. lb). At 
T — its imaginary part reads 

s ) - - I Im -lj + £do + £do (14) 
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It is seen that the first line of Eq. (14) represents mul- 
tiple Kondo resonances above Ef with the positions 
and weights given by Eqs. (11)-(13). The remainder of 
Eq. (14) describes the structure below Ef- The single- 
particle resonance (sdo) is strongly renormalized down- 
ward due to the logarithmic terms. The maxima of the 



log terms indicate weak peaks or often merely shoulders 
at the energies e' dm = —(idm — ido), i-e. at mirrorim- 
aged positions w.r.t. i d o- In contrast to the resonances 
above Ef, these mirrored peaks do not correspond to 
quasiparticle resonances, since there is no pole in Ada 
at these energies. Physically, they originate from spin 
fluctuations involving a virtual transition from an ex- 
cited state (m > 1) to the ground state (m — 0). This 
is only a virtual process, as the m > 1 states are not 
thermally occupied at T = [8]. With increasing T, 
< T < TKm, the peaks below Ef grow, as the m > 1 
states become thermally occupied. The multiple peak 
structure discussed above is shown in Fig. 2. Fig. 3 
shows fits of the NCA calculation of A d to photoemis- 
sion spectra of the Kondo system CeSi2 (M = 7) [12]. 
The experimental data show the appearance of low- 
energy peaks above Ef at elevated T, in quantitative 
agreement with the T dependence expected theoreti- 
cally from multiple Kondo resonances. 
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Fig. 3. Low-energy high resolution photoemission spectra of 
CcSi2- Open circles: Experiment; solid lines: NCA. Model pa- 
rameters determined from fits to the theory: Sdo— -1 -35 eV, SO 
splitting: 270 meV, CF splitting: 25 meV (lower quartet), 48 
meV (upper sextet); D — 3.7 eV; Tko ~ 35 K. The inset shows 
raw data before division by the Fermi-Dirac distribution [3]. 



3. Approximate conductance quantization 

We now show that in a symmetrical multi-level Qdot 
or a single-atom transistor (SAT) [13] at low T the 
linear response conductance G = d//dV(V = 0) is 
approximately quantized at 2e 2 /h [7], if the Qdot or 
SAT is in the Kondo regime, even though there may be 
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several transmission channels, m = 0, . . . , M — 1, and 
the lead-to-dot coupling matrix elements V m p are ran- 
dom. Such systems are described by the Hamiltonian 
Eq. (1), where, in addition, the conduction electron op- 
erators carry an index a = L,R indicating the left or 
right lead, and, hence the kinetic energy term reads, 
-Hkin = J2p<?J e P ~ M4^q<V<*, with fj, L/R = 0,V the 
chemical potentials in the left/right lead. 

The current through a left / right symmetrical device 
is [141 
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where Gda is the advanced local-orbital Green's func- 
tion and T the lead-dot couplings T^ nm , + r^ m , . G da 
and F are matrices in the space of local orbitals, 
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where E mm / (u) is the selfenergy due to intra-dot in- 
teractions. Since the ground state is a spin singlet [1], 
the quantum dot is at T <C Tk = min[TK m ] a pure 
potential scatterer for electrons traversing the system, 
and the following Fermi liquid relations hold [15], 
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The averaged electron number in the dot per spin, 
n,j,o-, can now be evaluated using the general relation 
^jln(G d _1 ) = (l-fj)-Gd and the Luttinger theorem 
Eq. (19), 
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This is a statement of the Friedel sum rule rida = 
i X/m ^mcr(O), since the scattering phase shift at the 
Fermi level in channel m is 8 ma (0) = arg[r-Gd CT (0)] mm . 
It may be re-expressed, using tr In Gda^ 1 = 
In det Gda-^ 1 , as 
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The scattering T-matrix of the device, T ■ Gda-, which 
determines the conductance G = dl /dV of the system 
via Eq. (15), is now evaluated by expressing the inverse 



matrix Eqs. (16), (17) in terms of its determinant, and, 
using the Fermi liquid property Eq. (18), we obtain at 
the Fermi energy (u> = 0, T < T K ) for M — 2, 



Im tr (r ■ Go (0)) = sin (7rn dCT ) + 
• ^ Re[det(ir-E'(0))] 



(21) 



ru(e d , 2 + £I, 2 (0)) + r22(£ d ,i + E'u(0)) ' 

If the transition amplitudes V mp are independent of the 
lead channels p, it follows directly from the definition 
ofT mm ' that the term oc sin(27rn [ j (T ) cancels. Eq. (21) is 
an exact result, valid for arbitrary microscopic param- 
eters T mm / , e dm , U , and n da . It is the generalization of 
the well-known unitarity rule of the single-level Ander- 
son impurity problem to the case of several impurity 
levels [7] . If there is at least one of the local levels sig- 
nificantly below the Fermi level (e d o < 0, |edo|/r mm / < 
1), a strong Coulomb repulsion U enforces n da ~ 1/2, 
implying via Eqs. (21), (15) a conductance close to the 
conductance unit, i.e. dl/dV fa 2e 2 /h (the factor 2 re- 
flects spin summation). 
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